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Abstract The influence of dynamical correlation effects on the magneto-optical properties 
of ferromagnetic Fe and Ni has been investigated. In addition the temperature 
dependence of the self-energy and its influence on the DOS and optical conduc- 
tivity is considered. Magneto-optical properties were calculated on the basis of 
the one-particle Green's function, which was obtained from the DMFT-SPTF 
procedure. It is shown that dynamical correlations play a rather important role 
in weakly correlated Fe and substantially change the spectra for moderately cor- 
related Ni. Magneto-optical properties obtained for both systems are found in 
better agreement with experiment than by conventional LDA calculations. 



1. Introduction 

Much information on the electronic structure of magnetic solids is gained 
by optical and magneto-optical measurements, being useful tools for analyz- 
ing the dispersion of (quasi-particle) bands. However, measured optical and 
magneto-optical spectra can hardly be interpreted without accompanying theo- 
retical calculations. For this purpose one in general has to solve a correspond- 
ing many-electron problem, which is impossible without the use of more or 
less severe approximations. For materials where the kinetic energy of the elec- 
trons is more important than the Coulomb interactions, the most successful first 
principles method is the Local (Spin-) Density Approximation (L(S)DA) to the 
Density Functional theory (DFT) [1], where the many -body problem is mapped 
onto a non-interacting system with a one-electron exchange-correlation poten- 
tial approximated by that of the homogeneous electron gas. For the last two 
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decades ab initio calculations of the optical and magneto-optical properties of 
solids based on thus approximation yielded a good basis for sucli an interpreta- 
tion, often leading to a quantitative agreement between theoretical and experi- 
mental spectra. The situation is very different when we consider more strongly 
correlated materials, (systems containing f and d electrons) since in all the 
calculations the LDA eigen-energies are implicitly interpreted to be the one- 
particle excitation energies of the system. It is well known that there are two 
possible sources of error connected with that approach: Firstly, the LDA pro- 
vides only an approximate expression for the (local) exchange-correlation po- 
tential. Secondly, even with the exact exchange-correlation potential at hand, 
one is left with the problem that there is no known correspondence between 
the Kohn-Sham eigen-energies and the one-particle excitation energies [2-5]. 

For an exact description of the excitation energies the non-local self-energy 
has to be considered. This, however, constitutes a many-body problem. There- 
fore, DFT-LDA calculations must be supplemented by many-body methods to 
arrive at a realistic description of the one-particle excitations in correlated sys- 
tems. To give an example, let us mention the GW approximation [6] which is 
well suited for the case of insulators and semi-conductors and has also been 
applied successfully to transition metals [6-9]. Another approach is to con- 
sider the Hubbard-type models where those Coulomb-interaction terms are 
included explicitly that are assumed to be treated insufficiently within DFT- 
LDA. Already the simplest Hartree-Fock like realization of such an approach 
called LDA-i-U [10] scheme allowed to improve considerably the description of 
the optical and magneto-optical spectra of strongly correlated systems (mostly 
containing rare earths elements [11, 12]). The main advantage of the LDA-hU 
scheme is the energy independence of the self-energy which allows to use 
only slightly modified standard band structure methods for calculating opti- 
cal and magneto-optical spectra. On the other hand the scheme works rather 
good only for extremely correlated systems, where Coulomb interactions (U) 
prevail considerably over the kinetic energy (bandwidth W). For moderately 
correlated systems (Uf«W) which applies for most 3d and 5/ elements and 
their compounds one has to take into account a non-Hermitian energy de- 
pendent self-energy to get a reasonable description of the electronic struc- 
ture. Nowadays there are several approaches available to deal with this sit- 
uation. The most advanced one is the Dynamical Mean-Field Theory (DMFT) 
[13]. DMFT is a successful approach to investigate strongly correlated sys- 
tems with local Coulomb interactions. It uses the band structure results calcu- 
lated, for example, within LDA approximation, as input and then missing elec- 
tronic correlations are introduced by mapping the lattice problem onto an ef- 
fective single-site problem which is equivalent to an Anderson impurity model 
[14]. Due to this equivalence a variety of approximative techniques have been 
used to solve the DMFT equations, such as Iterated Perturbation Theory (DPT) 
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[13, 15], Non-Crossing Approximation (NCA) [16, 17], numerical techniques 
like Quantum Monte Carlo simulations (QMC) [18], Exact Diagonalization 
(ED) [15, 19], Numerical Renormalization Group (NRG) [20], or Fluctuation 
Exchange (FLEX) [21-23]. The DMFT maps lattice models onto quantum 
impurity models subject to a self-consistency condition in such a way that the 
many-body problem for the crystal splits into a single-particle impurity prob- 
lem and a many-body problem of an effective atom. In fact, the DMFT, due 
to numerical and analytical techniques developed to solve the effective im- 
purity problem [13], is a very efficient and extensively used approximation 
for energy-dependent self energy S(i^). At present LDA-i-DMFT is the only 
available ah initio computational technique which is able to treat correlated 
electronic systems close to a Mott-Hubbard MIT (Metal-Insulator Transition), 
heavy fermions and /-electron systems. 

Concerning the calculation of the optical spectra we have to face the fol- 
lowing problem: one particle wave functions are not defined any more and the 
formaUsm has to applied in the Green function representation. Such a represen- 
tation has already been derived [24] and successfully applied for calculations 
in the framework of Korringa-Kohn-Rostoker (KKR) Green-function method 
for LSDA calculations. The only drawback of such an approach is that it is 
highly demanding as to both computational resources and computational time. 

In this paper we propose a simplified way to calculate optical and magneto- 
optical properties of solids in the Green function representation based on vari- 
ational methods of band structure calculations. 

The paper is organized as following: in section 2 the formalism for Green's 
function calculations of optical and magneto-optical properties that account 
for many-body effects through an effective self-energy is presented. Then, the 
DMFT-SPTF method for the calculation of the self-energy is considered. In 
section 3 the obtained results of our calculations for Fe and Ni are discussed 
and compared with experimental ones. The last section 4 contains the conclu- 
sion and an outlook. 

2. Green's function calculations of the conductivity tensor 

Optical properties of solids are conventionally described in terms of either 
the dielectric function or the optical conductivity tensor which are connected 
via the simple relationship: 

0-a/3(w) = -^(£a/3('^) " ^aji) ■ (1) 

The optical conductivity is connected directly to the other optical properties. 
For example, the Kerr rotation 6k{oj) and so-called Kerr elUpticity eK{oj) for 
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small angles and | Sxy Sxx \ can be calculated using the expression [25]: 



9k{oj) +i£K{oj) = 

(^xxil^) [l + ^(^xx{i^) 

The reflectivity coefficient r is given by 

_ (n - 1)^ + fc^ 
~ (n + 1)2 + A;2 



1/2 



(2) 



(3) 



with n and being the components of the complex refractive index, namely 
refractive and absorptive indices, respectively. They are coimected to the di- 
electric function via: 



n + ik = (exx + i£: 



xy) 



U/2 



(4) 



Microscopic calculations of the optical conductivity tensor are based on the 
Kubo Hnear response formalism [26]: 



1 



fiujV 



(5) 



involving the expectation value of the correlator of the electric current opera- 
tor J air). In the framework of the quasiparticle description of the excitation 
spectra of solids the formula can be rewritten in the spirit of the Greenwood 
approach and making use of the one-particle Green function G{E): 



(TapH = ^ / dE dE'fiE - - E') 

y J —oo J —oo 

Ti{j^QG{E')jpQGiE)} 



{E' -E + ir]){huj + E-E' + ir]) 

Tr {jpQGiE')j^QGiE)} 
{E' - E + ir]){huj + E' - E + irj) 



+ 



(6) 



where 9G(£') stands for the anti-Hermitian part of the Green's function, f{E) 
is the Fermi function and V is the volume of a sample. Taking the zero temper- 
ature limit and making use of the analytical properties of the Green's function 
one can get a simpler expression for the absorptive (anti-Hermitian) part of the 
conductivity tensor: 



1 i-Ef 
TTO; JEp-uj 



dEtr 



ja^G{E)jpQG{E + hu) 



(7) 
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The dispersive part of (Tapiyj) is connected to the absorptive one via a 
Kramers-Kronig relationship. 

The central quantity entering expression Eq.(7) is the one-particle Green's 
function defined as a solution of the equation: 

[HQ + t{E)-E]G{E) = i, (8) 

where is a one-particle Hamiltonian including the kinetic energy, the electron- 
ion Coulomb interaction and the Hartree potential, while the self-energy Ti{E) 
describes all static and dynamic effects of electron-electron exchange and cor- 
relations. The L(S)DA introduces the self-energy as a local, energy indepen- 
dent exchange-correlation potential Vxc{r)- As the introduction of such an 
additional potential does not change the properties of Hq we will incorporate 
this potential to Hlda and subtract this term from the self-energy operator. 
This means that the self energy S used in the following is meant to describe 
exchange and correlation effects not accounted for within LSDA. 

With a choice of the complete basis set the Green's function can be 
represented as: 

G{E) = Y.\i)GijiE){j\ , (9) 

ij 

with the Green's matrix Gij being defined as 

(10) 

Dealing with crystals one can make use of Bloch's theorem when choosing 
basic functions |ik)- This leads to the k-dependent Green's function matrix 

G^{E) = [Hf^ - EO% + Y1%{E)]-' . (11) 

Introducing the anti-Hermitian part of the Green's function matrix as 

g^iE) = '-[G%{E) - Gl{E)\ (12) 

and taking into account the above mentioned translational symmetry we obtain 
the following expression for the absorptive part of the optical conductivity: 

= — r I "^'^ E E)jfiiK E + hio) (13) 

TTLO J Ew—hiUi J 

with 



GiAE)= - E(:i\j) + (^\t{E)\j) 



Jt,{\^,E) = Y,QUE){n^\f\h) 

n 



(14) 
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The efficiency and accuracy of the approach is determined by the choice 
of |ik). One of the computationally most efficient variational methods is the 
Linear Muffin-Tin Orbitals method [27] which allows one to get a rather ac- 
curate description of the valence/conduction band in the range of about 1 Ry, 
which is enough for the calculations of the optical spectra (?iu! < 6 — 8 eV). 
This method has been used in the present work. A detailed description of the 
application of the above sketched approach in the framework of LMTO can be 
found elsewhere [28]. 

Calculation of the self-energy 

The key point for accounting of many-body correlations in the present ap- 
proach is the choice of approximation for the self-energy. As it was discussed 
in the Introduction one of the most elaborated modem approximation is DMFT. 

For the present work we have chosen one of the most computationally ef- 
ficient variants of DMFT: Spin polarized T-matrix plus fluctuation exchange 
(SPTF) approximation [23], which is based on the general many-body Hamil- 
tonian in the LDA-i-U scheme: 



where A = fm are the site number (i) and orbital (m) quantum numbers, 
a =T,i is the spin projection, c+,c are the Fermion creation and annihila- 
tion operators, Ht is the effective single-particle Hamiltonian from the LDA, 
corrected for the double-counting of average interactions among correlated 
electrons as it will be described below. The matrix elements of the screened 
Coulomb potential are defined in the standard way 



(12 34) = / drdr'^l{r)^^{r')v (r - r') V3(r)V'4(r'), (16) 



where we define for briefness Ai = 1 etc. A general SPTF scheme has been 
presented recently [23]. For d electrons in cubic structures where the one- 
site Green function is diagonal in orbital indices the general formalism can 
be simplified. First, the basic equation for the T-matrix which replaces the 
effective potential in the SPTF approach reads 



H 



Ht + Hu 



Ht 



= ^ J2 (^i^2|v|AiA^)c+^c+^,CA^^,CA;^, (15) 



(i3|t'^'^' (ifi)|24) = (13|i;|24) ^^^(13|v|56) x 




{iu) {in - iuj) (hQ ^T""' {in)\ 24^ , 



(17) 
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where uj = (2n + l)nT are the Matsubara frequencies for temperature T = 
(n = 0,±l,...). 

At first, we should take into account the "Hartree" and "Fock" diagrams 
with the replacement of the bare interaction by the T-matrix 

= ^EE(l3|^^"'(i^^)|23)Gr(^^^-^a;) 



a 3a' 



S£?M = -^^^ {13 \T-- mi 32) GUin-iu). 



Q 3 



(18) 



Now we rewrite the effective Hamiltonian (15) with the replacement (12 \v\ 34) 



by ^12 T'^'^ 3Aj in Hu. To consider the correlation effects described due to 

P-H channel we have to separate density (d) and magnetic (m) chaimels as in 
Ref.[21] 



di2 


















m^2 


= 4lC2V 





(19) 

Then the interaction Hamiltonian can be rewritten in the following matrix form 



(20) 



where * means the matrix multiplication with respect to the pairs of orbital 
indices, e.g. 



''22' 



34 



The supervector D is defined as 



D = (d, m°) , 



d+ 



and the effective interactions have the following form: 



12 

ll',22' 



2'V 
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<,22' = ^ E T^'' 1'2') - ^ E (12 2'1'> 



^n'?22' = oE^^'(l2 



l'2')-iE(l2|r""|2'l') 



''ll',22' — ''22', 11' — 



+ (12 



1'2') - (12 



1'2' 



12 



1'2' 



12 



1'2' 



2'l') + (12 



2'1' 



(21) 



To calculate the particle-hole (P-H) contribution to the electron self-energy 
we first have to write the expressions for the generalized susceptibilities, both 
transverse x"*" and longitudinal xH. One has 



[iuj) = [1 + * r^^ {iuo)\ * r^^ ( 



\IL0) 



(22) 



where 



-paa 
^ 12,34 



(r) = -G^(r)Gn-T)<5235i4 (23) 

is an "empty loop" susceptibility and T{iu)) is its Fourier transform, r is the 
imaginary time. The corresponding longitudinal susceptibility matrix has a 
more complicated form: 



and the matrix of the bare longitudinal susceptibiUty is 

x! 



1 / pTT _|_ pii pTT _ pii 



(24) 



(25) 



in the dd-, dmP-, mPd-, and mPmP- channels (d, m*^ = 1, 2 in the supermatrix 
indices). An important feature of these equations is the coupHng of longitu- 
dinal magnetic fluctuations and of density fluctuations. It is not present in the 
one-band Hubbard model due to the absence of the interaction of electrons with 
parallel spins. For this case Eqs. (22) and (24) coincide with the well-known 
result of Izuyama et. al. [29]. 

Now we can write the particle-hole contribution to the self-energy. Similar 
toRef.[22] one has 



(r) = E ^1M2 ir) Gi, (r) , 

34,0-' 



(26) 



with the P-H fluctuation potential matrix: 
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W^'^ {iuj) {iio) 



(27) 



were the spin-dependent effective potentials are defined as 




m 



rl- 
m 



(28) 



Here x")Xo dif^r from x",Xo replacement of T^^ ^ T^J^ in Eq.(25). 

We have subtracted the second-order contributions since they have already 
been taken into account in Eq.(18). 

Our final expression for the self energy is 



This formulation takes into account accurately spin-polaron effects because of 
the interaction with magnetic fluctuations [30, 31], the energy dependence of 
the T-matrix which is important for describing the satellite effects in Ni [32], 
contains exact second-order terms in v and is rigorous (because of the first 
term) for almost filled or almost empty bands. 

Since the LSDA Green's function already contains the average electron- 
electron interaction, in Eqs. (18) and (26) the static part of the self-energy 
S'^(O) is not included, i.e. we have 



3. Results and discussion 

The matrix elements of v appearing in Eq.(16) can be calculated in terms 
of two parameters - the averaged screened Coulomb interaction U and ex- 
change interaction J [23]. The screening of the exchange interaction is usually 
small and the value of J can be calculated directly. Moreover numeric cal- 
culations show that the value of J for all 3o? elements is practically the same 
and approximately equal to 0.9 eV. This value has been adopted for all our 
calculations presented here. At the same time direct Coulomb interaction un- 
dergoes substantional screening and one has to be extremely careful making 
the choice for this parameter. There are some prescription how one can get it 



(29) 



{iuj) = ^.''{iuj) - S^(0). 



(30) 
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Figure 1. The self-energy (a) of Fe for three different temperatures and corresponding densi- 
ties of states (b) and optical conductivities spectra (c). Full, dashed and dotted lines correspond 
to T = 125K, T = SOOK and T = 900iC, respectively. 



within constraint LDA calculation [2]. However, results obtained in this way 
depend noticeably on the choice of the basis functions, way of accounting for 
hybridization etc. Nevertheless the order of magnitude coming out from vari- 
ous approaches is the same giving the value of U in the range 1^ eV. In the 
present paper we are discussing the influence of the choice of U on the calcu- 
lated optical spectra. 

Another parameter entering SPTF equations is temperature. For a moment 
we are more interested in the low temperature properties while computationally 
the higher the tempreture is, less computationally demanding are the calcula- 
tions. This is why we decided first to consider the dependence of the self- 
energy on the temperature. 

In Fig. 1 we show the self-energy obtained for Fe for three different tem- 
peratures as well as corresponding densities of states and optical conductivities 
spectra. One can see that despite the differences in S are quite noticeable this 
leads only to moderate changes in the density of states and does not affect the 
optical conductivity. 

Much more important for the results is the parameter U. Fig. 2 shows as 
an example the real part t2g component of S for T = 300K in Fe for var- 
ious values of U. Despite the overall shape of the curve is practically the 
same the magnitude of the self-energy increasing with increase of U as it is 
expected from the analytical expressions. This change in self-energy leads to 
corresponding changes in the densities of states especially noticeable for the 
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Energy (eV) Energy (eV) 

Figure 2. The real part of t2g component of E for T — 300 AT in Fe for various values of U; 
left: spin-up, right: spin-down. 

minority spin subband. The influence of the choice of U on the optical proper- 
ties is even more pronounced (see Fig. 3). The low energy peak in the diagonal 
part of the optical conductivity shifts to the lower energies reaching the exper- 
imental position already for [7=1.5 eV. In the high energy part of the spectra 
large values of U lead to a structure around 5 eV not seen in experiment. 

Again, the value [7=1.5 eV gives also the best description for the shape of 
the experimental curve. (Note, that the experimental results for al.^{Lu) are 
multiplied by a factor of 1.7 to make the comparison more obvious.) A rather 
different situation occurs for the off-diagonal part of the optical conductivity. 




1 1 1 1 1 1 1 I 1 1 1 1 1 1 1 -0.6 1 1 1 1 1 H 

01234567 01234567 0123456 

Energy (eV) Energy (eV) Energy (eV) 



Figure 3. Optical conductivity (left: diagonal; middle: off-diagonal) and polar Kerr rotation 
(right) spectra in comparison with the experimental data of Fe. Experimental data for conduc- 
tivity are taken from Ref. [33, 34]; Kerr rotation spectra - from Ref. [35] 
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Figure 4. The real part of t2g self-energy for f/ = 1.5, 2, 3 eV and corresponding DOS plots 
for Ni; left: spin-up, right: spin-down. 



The low energy peak can be brought to the proper position only with [7=4 eV, 
at the same time the shape of the theoretical curves above 2.5 eV has a rather 
different structure in comparison with the experimental one, only crossing the 
zero axis at the same energy. However, a direct comparison of calculated a^y 
data with experimental ones may be somewhat misleading as experimentally 
this quantity cannot be measured directly and is usually obtained from ellip- 
sometric measurements and measurements of the Kerr rotation spectra. Thus 
in the left panel of the Fig. 3 we show our results for the calculated polar Ken- 
rotation spectra in comparison with experimental data. As one can see again 
the DMFT calculation with U = 1.5 eV describes the experimental data in a 
rather satisfactory way. 

If for Fe LSDA calculations already give a reasonable description of the op- 
tical properties and the many-body correlation effects, which improves only 
minor details, the situation in Ni is quite different. It is well-known that 
LSDA fails to describe the bandwidth for Ni, causing problems in the theoret- 
ical interpretation of all the spectroscopic experiments such as photoemission, 
x-ray emission, optics, etc. The main reason for this is the underestimation 
of electron-electron correlations which appear to be relatively strong in this 
metal. Again, as in the case of Fe, we carried out calculations with different 
values of U to find the best description of the spectral properties of Ni. In 
Fig. 4 we show the real part of the t2g self-energy for [/ = 1.5, 2, 3 eV as well 
as corresponding DOS plots. 

Despite the changes in the amplitude of the self-energy are important, all 
self-energies lead to rather small changes in the density of states, narrowing 
somewhat the bandwidth only and developing a low energy tail. Nevertheless 
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Figure 5. Optical conductivity and polar Kerr rotation spectra in comparison with the ex- 
perimental data of Ni. Experimental data for conductivity are taken from Ref. [36, 37]. Kerr 
rotation spectra - from Ref. [35] 



the diagonal part of the conductivity which reflects the convolution of the oc- 
cupied and unoccupied states is much more affected by the choice of U. The 
main change can be seen in the position of high energy peak which is placed 
by LDA about 1 eV higher in comparison with experiment. Accounting for the 
correlation effects shifts this maximum bringing it to the proper position for 
U = 3 eV. The low energy part of the spectra does not reflect too much influ- 
ence of the U parameter and deviates just slightly from the experimental curve. 
For the off-diagonal part of the conductivity an improvement as compared to 
LDA is not so pronounced as for the diagonal one, though the spectra getting 
closer to experiment. It is worth to note that the actual value of U doesn't 
change the calculated spectrum of ioa'^y{io). But again, as mentioned in the 
case of Fe it is worth to compare calculations with directly measured Kerr ro- 
tation spectra presented in Fig. 5. As one can see, the improvement compared 
to LSDA results is substantional but our results are still far from experiment 
concerning the peak position both in the infrared and visible parts of the spec- 
tra. This disagreement is apparently coming from the approximation that has 
been made and is much more pronounced in the off-diagonal part of conduc- 
tivity as it is more sensitive to the details of the electronic structure being the 
result of complex interplay of exchange splitting and spin-orbit coupling. 

It is still unclear whether the mentioned problems are coming from the 
single-site approximation for the self-energy (DMFT) itself or whether they 
are reflecting the limitations of the simplified FLEX method of solving the im- 
purity many-body problem. To find out an answer more elaborated solvers like 
QMC have to be used. 
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4. Conclusion and outlook 

In the present paper we show a way to account for the particle-particle corre- 
lations in the theoretical description of optical and magneto-optical properties 
of the ferromagnetic 3d metals. We show that the dynamical correlations play 
an important role even in weakly correlated materials like Fe and can sub- 
stantially change the shape of the spectra for moderately correlated Ni. Even 
a rather simple way of accounting for dynamic correlation allows to improve 
theoretical results substantionally though not giving the perfect agreement with 
experiment. 

Thus to go further one has to use more elaborated technique to obtain the 
self-energy both within DMFT and beyond (for example, new DMFT-i-GW 
approximation). Work along this line is in progress. 
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